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Abstract 



We show that the Hamiltonian of four-dimensional Lorentzian gravity, defined on a space 
of real, SU(2)-valued connections, in spite of its non-polynomiality possesses a natural quan- 
tum analogue in a lattice-discretized formulation of the theory. This opens the way for a 
systematic investigation of its spectrum. The unambiguous and well-defined scalar product 
is that of the SU(2)-gauge theory. We also comment on various aspects of the continuum 
theory. 
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In spite of considerable advances in our understanding of the canonical quantization of 
gravity, based on the use of complex connection variables [1] and quantum representations 
on spaces of Wilson loops [2], some basic problems have remained unsolved. Firstly, despite 
the apparent simplicity and polynomiality of the Hamiltonian in this approach, even at a 
formal level only a few (and for the most part physically uninteresting) solutions to the 
Wheeler-DeWitt equation have been found in any loop representation. Secondly, treatment 
of the reality conditions (which have to be imposed on the SL(2, (D)-valued basic variables to 
recover real Einsteinian gravity) continues to be troublesome in the quantum theory. (The 
heat kernel measure on SL(2, C)-wave functions used in [3,4] provides a kinematic scalar 
product for complex Wilson loops (modulo a possible hidden metric dependence), but does 
not incorporate correctly the reality conditions for Ashtekar gravity.) 

Faced with these difficulties, it may be time to remember that there exists a version of 
Hamiltonian gravity in terms of real connection variables [5] (for a corresponding action prin- 
ciple, see [6]), and to re-evaluate the achievements and drawbacks of the complex formulation 
as compared to this real alternative. There is a "unified" derivation of the two connection 
representations in the classical theory: starting from the 3+1 formulation in terms of a cano- 
nically conjugate pair (P", K l a ) of S'0(3)-valued variables, where Pf denotes a dreibein (with 
density weight 1) and K % a = -j=K a iP bl is (whenever the Gauss law constraints are satisfied) 
the extrinsic curvature with one index raised, one may define a canonical transformation 

E- = olPZ 

■ (1) 

Ai = r a + pK, 

where T = T(P) is the spin connection compatible with P, and a and (3 are two non- vanishing 
constants. In terms of the new variables, one has 

{Ai(x), E b (y)} = a(3{Ki(x), P b {y)} = -af55)5 b j\x - y). 
The Hamiltonian constraint function in terms of (A, E) reads 

e i jkE a EbjF k b _ { ?_ + 2)E \a E b] _ ^ ){Ai _ (2) 

Ashtekar's choice is a = 1, (3 = —i, which makes the second, non-polynomial term drop out of 
the Hamiltonian. The drawback is that according to formula (1), the connection variable A\ 
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is now complex. Our real version of connection gravity consists in choosing a = 1 and (3 = —1 
[5]. The functional forms for the remaining spatial diffeomorphism constraints E^F^ b = 
and Gauss law constraints V a £"f = are independent of the choice of a and (5. (As usual, 
we denote the field strength of the connection A\ by F and its covariant derivative by V.) 
Their treatment therefore does not have to be changed with respect to the usual, complex 
formulation. In particular, we will continue to use Wilson loop variables in the quantum 
theory. In the real formulation, (2) becomes 

retaining a "potential" term (which is a misnomer since it depends both on coordinates and 
momenta). Starting from a form equivalent to (2) (formula (14) in Barbero's paper [5]) for 
H pot , one finds after some algebra that it may be re-expressed as a polynomial in A and E, 
up to determinantal factors, namely, 

(4) 

and up to terms proportional to the Gauss law constraints. 

Let us now recall some features of the complex formulation with the Hamiltonian = 
e l i k EfEjF^ b . There exist simple solutions to H®ty = in the loop representation, both in 
the formal continuum approach [2] and on the lattice [4], where the loop state \I> depends on 
(smooth) non-intersecting loops (see also [7] for a generalization within a lattice language). 
This is a straightforward consequence of the antisymmetry of H® in the spatial indices a 
and b. However, these solutions are probably not interesting from a physical point of view, 
because they correspond to zero-eigenstates of the volume operator [8,9]. We are not aware 
of any non-trivial (in this sense) solutions that have been found by tackling the equation 
H^qr = directly. 

Other interesting features of the complex formulation are that all four diffeomorphism 
constraints can be solved by making a so-called Capovilla-Dell-Jacobson ansatz [10], and 
that there exists a formal solution to the Wheeler-DeWitt equation (with a specific factor 
ordering for H ), given by the exponential of the Chern-Simons action for the complex 
connection A [11]. However, in the absence of a proper treatment of the reality conditions, 
the significance of these properties for the full gravitational theory has remained unclear. 
There is a version of the CDJ-ansatz in the real theory [12], but the remaining Gauss law 
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constraints are considerably more difficult. The real Chern-Simons term does not seem to 
play a special role (apart from being a generating functional for the -B-field, B at = ^x] ahc F^ 
like in the complex theory). This also implies that the solutions to the Hamiltonian related 
to knot invariants (obtained by a formal loop transform of the Chern-Simons term from the 
connection to the loop representation [13]) do not carry over to the real theory. The absence 
of these "nice" features from the real theory may lead one to wonder whether they may not 
go away also in the complex case once reality conditions are properly taken into account. 

The (quantum) Hamiltonian of the complex theory has a natural representation in 
terms of loop functions, since the components of the field strength F both in the continuum 
[2] and on the lattice [16] can be obtained by considering infinitesimal planar Wilson loops. 
We will show below that a similar statement holds for the real Hamiltonian H 1 ^, i.e. in spite 
of its non-polynomiality, the potential term has a natural representation on quantum loop 
states, at least in the lattice formulation. This makes the search for zero eigenvectors of 
accessible numerically. Clearly the potential term presents a computational complication, 
which however in the absence of any explicit solutions of H^ty = seems to be a matter of 
degree rather than one of principle. The big advantage of the real formulation is the presence 
of a well-defined and unambiguous scalar product on Wilson loop functional in the quantum 
theory, that induced by the Haar measure of SU(2). 

The following discussion will take place within the discrete lattice framework, with oc- 
casional comments on a possible continuum formulation. Furthermore, we will focus on the 
discussion of the potential term, which is new with respect to previous treatments [16,4,7]. 
It may be worthwhile noticing that in the lattice approach, both the state space and the 
operators get regulated simultaneously, since they share the same support (in terms of lattice 
links); in discretizing the state space, we get a regularization of the Hamiltonian "for free". 

Recall now the basic ingredients of the Hamiltonian lattice formulation for theories based 
on a space of connections [17]. Our lattice will be a cubic N x N x TV-lattice, with periodic 
boundary conditions. The basic operators associated with each lattice link I are in our case 
an SU (2)-\mk holonomy V (represented by multiplication by V), together with its inverse 
V~ l , and a pair of canonical momentum operators pf and p~ , where i is an adjoint index. 
The operator pf(n, a) is based at the vertex n, and is associated with the link I oriented in 
the positive d-direction. By contrast, p~(n + la, a) is based at the vertex displaced by one 
lattice unit in the d-direction, and associated with the inverse link l~ 1 (a) = /(— d). The wave 
functions are elements of XiL 2 (SU(2),dg), with the product taken over all links, and the 
Haar measure dg. The basic commutators are 
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[V A B (n,a),V c v (m,b)] =0 



[p+(n,a),V A c (m,b)} = - l -5 nm 5^T iA B V B C 

[p-(n, a), V A C (m, b)} = - % - 5 nm 5^ V A B r iB c ( 5 ) 
\pt{n,a),pf(m,b)] = ±i 5 nm 6^ b e ijk p k 
\pf(n,a),pj(m,b)] = 0, 

where e^ k are the structure constants of 577(2). Our normalization for the 5C/(2)-generators 
Tj is such that [Ti,Tj] = eij k T k and A a = A l a Ti/2. Taking into account the expansions 



V A B {b) = t A B + aA bA B + 0{a 2 ) 
pf(b) = a 2 E b l+ 0(a 3 ) 

of the corresponding classical quantities for small lattice spacing a, one derives the following 
expansion 



-i'&(r i K(n,S)p^(n + i 8 ,c)r J -F(n,S)- 1 )-pi t (n,c) = 

a 3 ^? + e^Ay^) + 0(a 4 ) = a 3 V b E? + 0(a 4 ). 

The prefactor —1/2 occurs because of Trr^Tj = — 2<5^ . 

For reasons of symmetry we will from now on use the averaged momenta pi := (pf +p~) /2 
and their quantum versions. Motivated by (7), we may represent the quantum covariant 
derivative V b Ef by the lattice operator 



--Tr ( n V(n, b) Pj (n + i g , c) Tj V(n, 6)" 1 ) - fc(n, c). (8) 



Note, however, that this operator is well defined only on those Wilson loop states that for 
each occupied link l(n,c) have also the neighbouring "parallel" link l(n + l b ,c) occupied. 
(A lattice Wilson loop is a gauge-invariant function of the form Tr V(li)V(l2)---V(l k ), with 
7 = li o l 2 o ... o l k a closed loop of lattice links.) 

This happens because the left-hand side of expression (7) is a finite difference. If one of 
pj(n + lg, c) and pi(n, c) but not the other vanishes on a state the result of the action of 



(8) on ^ for small lattice spacing a is of lower order in a and diverges in the limit as a — ► 0. 
In particular, loop states with only sparse intersections have a singular behaviour under the 
action of the covariant derivative operator (8). Moreover, since we do not want to distinguish 
any particular direction on the lattice, we will instead of (8) use the averaged version 

- J(Tr (nV(n, S)p,-(n + ig, c) Tj V(n, 

- Tr { Tl V{n - i g) br'pjin - ig, c) Tj V(n - ig, 6))). 

A remark similar to the one made above concerning the allowed loop states applies to this 
operator as well. 

To obtain a well-defined lattice operator H pot , we still must take care about the de- 
terminantal factor (det.E) -2 . As shown in [9,14], det.E = ^ i r] a f }C e ljk EfEjE^ possesses a 
natural quantum lattice analogue D(n) := ^7] a b c e l i k pi(n,a)pj(n,b)pk(n,c). Since the latter 
is a selfadjoint operator, there exists a Hilbert space basis in which it is diagonal. For the 
gauge-invariant sector, this basis is most easily constructed in terms of so-called spin network 
states, certain (anti-)symmetrized, real linear combinations of Wilson loop states. 

A spin network associates a positive "occupation number" with each lattice link, which 
may be interpreted as counting the number of (unoriented) flux lines of basic spin-^ repre- 
sentations along the link, and also keeps track of the way in which those flux lines can be 
contracted gauge- invariantly at the vertices (see [18] for more details). A concrete way of 
constructing elements of the spin network basis is to begin with sets of Wilson loops with 
fixed occupation numbers and arbitrary intertwiners, and then select linearly independent 
sets of intertwiners at the vertices (which generally, in terms of a loop language, are still 
related by Mandelstam constraints). 

In terms of such states, the diagonalization of the operators D(n) is reduced to the 
diagonalization within finite-dimensional subspaces of the Hilbert space. In the resulting 
diagonal basis we can meaningfully define quantum representations of arbitrary functions of 
(det E) in terms of their eigenvalues. (Investigations of the spectrum of (two related but not 
identical versions of) the volume operator have been performed in [14,15].) In particular, if we 
restrict ourselves to eigenstates with non- vanishing eigenvalues, we can quantize (det E)~ 2 on 
the lattice. There is no immediate analogue of this construction in the continuum, although 
one can define a quantized version of the classical volume function J ^\ det E\ [8] . Since the 
naive local quantum operator det E vanishes at all points of a loop state without intersections, 
1/det E is ill-defined almost everywhere on a typical loop state. In the continuum there thus 
seem to be no good analogues of the lattice states with "volume everywhere". A way out 
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of this may be to either use smeared-out wave functions and/or regularize the quantum 
operators appropriately. Alternatively, one may multiply the Hamiltonian by a factor 
(det E) 2 . Classically, this changes the constraint algebra at most by terms proportional to the 
constraints. In particular, the commutator of two Hamiltonians is just rescaled by a factor 
of (det.E) 4 . Quantum-mechanically, however, operator ordering problems may appear. 

Possibly the problem is not as serious as it seems since along smooth pieces of loop where 
det E vanishes, also (det E) 2 H pot does, and one may be able to regularize the action of H pot 
to a finite value. Another problem, also relevant to the lattice approach, is that detE has 
many zero- volume eigenstates even at loop intersections [9,14]. Hence the question is whether 
one can consistently restrict the Hilbert space so that the action of H pot is always well defined 
(recall that [H^,detE] ^ 0). For example, no immediate problems arise if one chooses a 
factor ordering for H pot with the det E-terms to the right and discards all zero- volume states 
by hand. 

To summarize: there is a well-defined regularized, self-adjoint operator expression for 
the real Hamiltonian (3) on the lattice, at least on a large subsector of the Hilbert space of 
gauge-invariant functions. This operator is subject to the usual ambiguities with regard to 
factor ordering and addition of higher-order terms in the lattice spacing a. For example, we 
may choose 

H n ' latt (n) = -e ijk Tr (V(n,P ai )T k ) Pi (n, a) Pj (n, b) 

Vbib 3 b A {Vk{n, a 3 )pi(n, a 4 )p m (n, b 3 )p n (n, 6 4 ) 

- 2p m (n,a 3 )p n (n,a A )p k (n,b 3 )pi(n,b A ) ) p k (n, a 2 )p m (n, b 2 ) x 

jTi(TiV(n - l a2 ,a 2 )~ 1 p s (n - la 2 , ai)T s V{n - la 2 ,a 2 ) (10) 

- TiV(n,a 2 )p s (n + la 2 , ai)T s V(n, a 2 ) _1 ) x 
Tv(T n V(n- l^Mr^tin- i k M)T t V(n - i &2 ,6 2 ) 

- T n V(n, b 2 )p t {n + l h M)T t V{n, h)- 1 ) D{n)- 2 

for the Hamiltonian localized around a vertex n. In (10), all spatial indices a, b etc. are 
summed over, and V(n, P^) is the holonomy associated with a plaquette loop in the a-b- 
plane. 

Note that the Hamiltonian of metric gravity (also containing inverse powers of det g = 
\detE\) cannot be treated in a similar way. The construction above depended on (i) the 
reformulation of canonical gravity in terms of connection variables, hence (ii) the possibility of 
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choosing a gauge-invariant Hilbert space of Wilson loops, therefore (iii) the diagonalization of 
the operators D(n) in terms of spin network states, together with (iv) a natural regularization 
of the covariant derivative terms in H pot . 

Having thus set the stage for a systematic investigation of the Hamiltonian eigenvalue 
problem H*> latt y = 0, we will now describe some technical problems that have to be ad- 
dressed for its solution. Consider the action of a local lattice Hamiltonian H^ ,latt {n) on a 
spin network state ^. Since the momentum operators do not change the occupation num- 
bers ji of links, it a priori looks as if this action would result in a set of loop states with 
Aji £ {0, +1, +2}, depending on the contributions V(link) to the various links in a neighbour- 
hood of n, coming from the kinetic and potential parts of the Hamiltonian. Unfortunately, 
life is not as simple. 

Take, for example, the action of the polynomial part O of H pot on a spin network state. 
Our computations show that the resulting terms generically do not form a set of states that 
combine in a simple way to give one (or a small number of) spin networks, because the 
operator action does not preserve the total (anti-)symmetry over link permutations of the 
spin network. The result must always be expressible as a unique linear combination of spin 
network states, but it turns out that this decomposition in general contains states whose 
occupation numbers differ from those of the original state within a whole range of values. 
For instance, the action of O(n) on a state ^ with occupation number ji for some link 1^ 
based at the vertex n may result in a sum of spin network states with jVs taking any one 
of the positive values ji + 4,ji + 2,ji,ji — 2, — Moreover, through retracings of the form 
V(n' , a)V(n' , a) -1 = 1 occurring during the decomposition of into independent spin 
networks, even links may be affected that where not acted on directly by 0(n) in the first 
place. This reveals a somewhat unpleasant property of the spin network states which in 
a sense are "maximally non-local" (as opposed to sets of maximally localized Wilson loop 
functions that one may favour in certain gauge-theoretic applications [19]), especially in 
conjunction with our requirement of selecting only spin networks with non-zero volume at 
every vertex. 

We therefore conclude that the investigation of the spectrum of H^ latt requires the pre- 
sence of an efficient algorithm for generating independent spin network states and computing 
inner products of such states. We reckon that even in the discretized lattice version the 
spectral problem is sufficiently complicated so as to make further approximations necessary. 
Since we can calculate matrix elements of the Hamiltonian explicitly, we can neglect small 
contributions, depending on suitable perturbation parameters such as those characterizing 
the spin network states or related to the bare gravitational coupling constant G. More details 
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on our investigation of these issues will appear elsewhere. 

Coming back to the continuum theory, one can show that the linearized limit of the real 
connection formulation coincides with that of the complex one (c.f. the treatment of usual 
Ashtekar gravity in [20]). In the real case, the linearized versions of the kinetic and potential 
terms in (3) become proportional to each other and add up to the expected result. Likewise, 
the large-G limit, as, for example, discussed in [21], is unchanged. For the G — > limit [22], 
this does not necessarily seem to be the case. 

Our real treatment (or an appropriate continuous analogue) can be viewed as comple- 
mentary to another approach that has recently been suggested for dealing with the complex 
version of the theory, namely, making use of a "generalized Wick transform" [23]. In this 
ansatz, one tries to define a transformation W = expC, with C = ^ f K^Ef, between two 
Hilbert spaces where in one the Hamiltonian has the simple form H E = e i:jk EfEjF^ b and in 
the other the more complicated form of the real theory. As far as we understand, the diffi- 
culties in making the operator W well defined in the continuum quantum theory are roughly 
comparable to those of constructing the continuum Hamiltonian H^. (Note that, like the 
phase space functional C, also the potential term H pot can be written in terms of Poisson 
commutators of the quantities J(det E)~ 2 H E , f \J det E, A\ and Ef.) However, even if 
these could be overcome, the problem of finding non-trivial solutions to the Wheeler-DeWitt 
equation would still remain. 

Acknowledgement. The author is indebted to F. Barbero for numerous discussions on real 
connection gravity. 
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